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Abstract 

We consider the effects of a quench to T = on a spin system with axial next-next 
nearest neighbour Ising interactions, evolving under a conserved 3-spin flip dynamics. 
Such a model is motivated by the kinetics of stacking layers in polytypes near the 
3C — 6H transition. We find that the system generically gets arrested in interesting 
metastable states which have inhomogeneously distributed quiescent and active regions. 
In such arrested states, the autocorrelation function decays as a stretched exponential 
~ exp{—{t/T )z). The latter feature can be understood in terms of a mapping of the 
dynamics within active stretches to the well known simple exclusion process of particles 
on a line, and bounds can be put on r . 
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1 Introduction 

When a statistical system is quenched by suddenly reducing the temperature T, we encounter 
an interesting situation, namely the nonequilibrium evolution of a high-T state in a low-T 
environment. Two outcomes are possible. The system may approach the equilibrium ground 
state, perhaps through a coarsening route. Alternatively, it may get arrested in a long-lived 
nonequilibrium state whose characteristics depend both on the interactions in the system 
and on gross features of the dynamics, such as conservation laws. In this paper, we examine 
the consequences of quenching a 1-d Ising spin system with competing interactions, which 
evolves by a multispin-flip dynamics which conserves an infinite number of quantities. Both 
the interactions and the kinetics are motivated by the use of the extended axial next nearest 
neighbour Ising (ANNNI) model to describe polytypism in close-packed crystals, as explained 
below. 

Within this model, we find that the nature of the arrested state obtained by a quench to 
T = depends on the degree of competition in the interactions. In several cases, the state 
is one in which there are dynamically active and quiescent regions interspersed randomly 
through the system. Such inhomogeneously active and quiescent (IQA) states have inter- 
esting dynamical properties, and are the focus of this paper. In particular, we show that in 
this state, the autocorrelation function follows a stretched exponential decay law, and that 
this can be related to the dynamics of the simple exclusion process in finite regions. 

2 Polytype Kinetics and the ANNNI Model 

Materials such as silicon carbide and zinc sulphide are known to undergo a number of tran- 
sitions from one close-packed arrangement to another when external conditions such as the 
temperature and pressure are changed |I|]. Each crystalline structure is associated with a 
particular sequence of stacking close-packed triangular layers; transitions between different 
polytypes correspond to rearrangements of the stacking sequence. 

Let us denote the three possible arrangements of the triangular lattice in a layer by A, B 
and C ||. Stacking of layers are then described by sequences such as ABCABC- • • (3C, face- 
centered cubic), ABCACB- • ■ (6H, hexagonal close-packed) and ABAB- • • (2H, hexagonal 
close-packed), each a different polytype. The polytypic constraint of close-packing implies 
that no two successive layers can be of the same type. 

The extended ANNNI model with third neighbour interactions plays an important role 
in understanding the occurrence of polytypes 0, [| . The Hamiltonian is 

Ti- = —J\ 2^ S{Si + \ — J2 2_^ Si.Si + 2 — J3 2_^ SiSi + 3 (1) 

i i i 

where Ji, J 2 and J 3 are coupling constants, competition between which can lead to the 
formation of long-period structures. The Ising spins Si = ±1 in Eq.(l) are related to 
successive pairs in the stacking sequence of A's, B's and C's as follows: Pairs AB, BC, CA 
correspond to Si = 1, while pairs BA, CB, AC correspond to Si = — 1. Thus for instance 
the 3C structure is translated to the spin sequence X T T ' ' '> the 6H structure is TTTIII 
while the 2H structure is tltl • • •• A portion of the ground state phase diagram is shown in 
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Figure 1: A portion of the ground state phase diagram of the third neighbour ANNNI 
model. The lines are first order boundaries; thick lines denote multiphase boundaries with 
an exponentially large number of degenerate ground states. 



Fig. 1. The phase boundaries shown all correspond to first order transitions; the heavy lines 
denote boundaries on which an exponentially large number of ground state configurations 
are degenerate ||. 

Besides providing a description of polytype phases, the extended ANNNI model can also 
be used to model polytype kinetics. The first such study was concerned with the kinetic 
effects of quenching across the T = 2H-6H phase boundary ||, and the nature of the 
resulting arrested states. Here we study the effects of quenching the model from T = oo to 
T = in the vicinity of the 3C-6H phase boundary. While modeling kinetics, it is important 
to allow for only those spin re- arrangements which correspond to physically admissible moves 
of the stacked layers. While local changes in the stacking sequence correspond to local 
changes in the spin configuration, the reverse is not necessarily true. For instance, a single 
spin flip corresponds to a deformation fault which involves the simultaneous movement of a 
macroscopic number of layers. 

Near the 3C-6H transition, there is experimental evidence |6|, |7| that the predominant 
kinetic move is a double-layer displacement which involves the simultaneous movement of 
two adjacent layers, for example, ■ ■ • CABC ■■■—>■■■■ CBAC • ■ ■ ; two successive layers in the 
sequence can interchange provided that such an interchange does not violate the polytypic 
constraint of no successive occurrences of the same letter. In spin language, the double-layer- 
displacement move corresponds to flipping triplets of successive spins, e.g. ■•• TT-UJ-TI ~* 
TTTTTTI ~~ * TIIITTT • • •• This is a generalization of the familiar single-spin-flip Glauber 
dynamics, but there is an important qualitative difference between the two. While Glauber 
dynamics in nonconserving, it turns out that triplet flip dynamics leaves invariant an infinite 
number of conserved quantities [§, |9|]. These conservation laws have important consequences 
for time-dependent properties. 
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Figure 2: Schematic depiction of partitioning of phase space under trimer-flip dynamics at 
T = oo and T = 0. Configurations which can be reached from each other through the 
dynamics are shown with the same symbol. The set of dynamically connected configurations 
defines a sector. The additional constraints that operate at T = lead to a reduction in the 
number of configurations as well as the number of sectors. 

3 Quenching from T = oo to T = : the IQA State 

Our problem involves studying the effect of quenching to T = a system with competing 
interactions and conservation laws. The dynamics consists of stochastic flips of triplets, and 
the T = environment allows only moves with energy AE < 0. 

Figure 2 is a schematic depiction of the available phase space (i) at T = oo, and (ii) along 
the 3C — 6H multiphase phase boundary at T = 0, where any configuration with parallel 
spin stretches with 3 or more spins is an allowed ground state. At T = oo, it is known 
||. that the infinity of conservation laws leads to a partitioning of the 2 L -large phase space 
into a large number Ni of sectors; N^ grows as ~ fi L where /i ~ 1.618 is the golden mean. 
In case (ii), by contrast, using a transfer matrix method we find that the total number of 
configurations is reduced to ~ k l with k ~ 1.466. This ground state subspace is further 
partitioned into Nl ~ j^ L sectors, with \x ~ 1.167 as a consequence of the infinitely many 
conservation laws. 

Evolving the system from a random initial configuration, with the AE < condition 
enforced, corresponds to a rapid quench from T = oo to T = 0. The energy non-raising 
condition imposes local constraints on whether or not a pair of chosen spins can actually be 
flipped; these constraints are a function of J±, J2 and J3. The normalized energy changes 
Ae = AE/Ji involved in flipping a triplet | j ] to ||| depend on the environments of the 
triplet, and are given in Table 1. There are 20 distinct local environments; the other unlisted 
environments (out of a total of 64 possibilities) are related to these by reflection symmetries. 



Table 1: Energy changes for triplet flips in various environments 



Dynamica 


moves 


Ae 


trttTTT-flDt -> 


ftWUftftft 


4(1 + 2j 2 + 3j 3 ) 


Ir-lr-lrTTT-U-U-U -> 


1HKUUW 





tftfrTTTtfN)- -> 


WIUW 


4(l + 2j 2 + 2j 3 ) 


tfKrTTTWT - 


fttfrJ-UW 


4(l+J 2 + 2j 3 ) 


WrTTTW -> 


Mtmw 


4(1+J2+J 3 ) 


1r-fr1rTTT^1rir -> 


fttfr-iuwr 


4(j 2 + 2j 3 ) 


MtTTTW -> 


wmw 


4(j2+J 3 ) 


IKKrTTWt -> 


wmw 


4j3 


1rMTTT1rM -> 


wmtM 


4(j2+J 3 ) 


1rMTTT1NMr -> 


wwwr 


4j 3 


ftWTTW -> 


www 





1rMTTTM1r -> 


WWW 


4(-l+j 3 ) 


WTTTW -> 


WWW 


-4 


IrMTTTWr -> 


WWW 


-4(1+ j 2 ) 


WrTTTW -> 


WWW 


4(1+J 2 +J 3 ) 


WrTtmt -> 


www 


4(1 + h) 


MtTTTW -> 


www 





WrTTWt -> 


WWW 


-4j 2 


1rWTT1rM -> 


www 





WTTWft - 


WWW 


-4(l + 2j2+J 3 ) 



For all values of j 2 = Jij J\ and j 3 = J3/J1 for which the Ae's listed in Table 1 have the 
same sign, the system will reach an arrested state of the same character. 

We used numerical simulations to investigate the properties of arrested states obtained by 
quenching to points on the 3C — 6H boundary and in the 3C phase (points P and Q in Fig. 
1, respectively). In both cases, the quench does not lead to the ground state. For instance, 
at point P, the system does not approach the subspace depicted in Fig. 2, but rather an 
arrested state which is infinitely long-lived at T = 0. Figure 3 gives an idea of the nature 
of this arrested state. The state has alternating quiescent and active stretches of random 
lengths. The lengths of active stretches have an exponential distribution, Ap(£) ~ exp(—Xp£) 
with Xp ~ 0.22. We monitored the decay of the spin autocorrelation function Cp(t) in this 
state (Fig. 4), and found that it decays as a stretched exponential ~ exp(— (£/rp) 1//3 ). We 
repeated the quench for a point Q (Fig. 1) in the 3C phase, away from the multiphase line. 
The equilibrium ground state is j ] | j • • • in the case Q, in contrast to the large degeneracy of 
ground states at P. Nevertheless, on quenching towards Q, we found again an IQA state with 
an exponentially distributed number of active stretches, Aq(£) ~ exp(— Xq£) with Xq ~ 0.13. 
Further, the autocorrelation function Cq(£) decays as ~ exp(— (t/rq) 1 / 3 ), but with Tq > Tp 
(Fig. 4). A study of the nature of the excitations in the active stretches in the two cases P 
and Q shows that they are quite different in character. Despite this, as argued below, they 
show the common feature of diffusing excitations confined to randomly distributed stretches 
of finite lengths, and this leads to the same form of stretched exponential decay of C(t) in 
both cases. 
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Figure 3: The space-time diagram of an IQA arrested state (left) with activity depicted in 
white. The activity pattern in a single active stretch of the IQA state (right) shows the 
diffusive motion of active walkers, with hard core interactions. 

Quench to P (j'2 = 0, j'3 = —1/3, on 3C — 6H boundary): In this case, the active 
regions consist of parallel spin patches (PSPs) in a background which is one of two types: 
(a) antiparallel doublets TTIITT ' ' ' (b) antiparallel doublets with some || pairs interspersed 
between doublets, e.g. TTITIITT ' ' '• Flipping a triplet at the edge of a PSP is an allowed 
(AE = 0) move, and results in a translation of the PSP by either one or two lattice spacings, 
e -g- TTIITTTTTII - ►TTiiiiiTTii- m the process, the size of the PSP may fluctuate between 
4 and 5, depending on the surroundings, e.g. TTITIIIIITT^TTIT 1 TTIITT- Backgrounds of 
type (b) can support only a single PSP, as two PSPs around an Tl pair lead to a frozen 
structure. By contrast, backgrounds of type (a) can support several mobile PSPs. In this 
case, two PSPs have a hard core repulsion and cannot pass through each other. These 
statements can be easily checked using Table 1, since AE > moves are disallowed. Thus 
we have a gas of extended objects with hard core interactions, diffusing in a confined region. 
The boundaries of quiescent regions have configurations only of a few types — for example, 
one such left-boundary is TTIII •••> which is stable with respect to the adjacent active 
stretches described above (again specific moves can be checked using Table 1). 

Quench to Q (j 2 = 1, js = 0, in 3C phase): In this case, active regions consist of stretches 
of ground state configurations, separated by domain walls, e.g. ■ ■ • TTTT I lllll I TTTT ' ' '• 
The spacing between successive domain walls is (3m + 2) with m = 0,1,2,.... A triplet 
adjacent to a wall can flip, whereupon the wall moves by three lattice spacings. Two diffusing 
walls exhibit a short range repulsive interaction, as they can approach no closer than two 
lattice spacings. Once again, the dynamics is tantamount to a gas of objects with hard core 
repulsion diffusing in a confined region. The quiescent boundaries in this case are also of a 
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Figure 4: Autocorrelation functions C(t) for the two types of IQA arrested states at the two 
points Q (filled circles) and P (empty triangles) decay as exp(—(t/T*)3). 



few types, e.g. a stable left-boundary of a quiescent region is TTiTT •••> it can De checked 
that any attempt to flip its edge spins leads to energy raising. 

In the active stretches in both cases P and Q, we have diffusing walkers with hard core 
repulsion confined to finite regions between impenetrable walls. Let the number density of 
such walkers in any active stretch be p. The autocorrelation function Ci (p, t) (where X 
stands for either P or Q), of such a stretch is expected to decay at large times as exp(—t/r x ). 
Here Tx is the relaxation time which depends on the size I of the finite stretch. A crucial 
point to note is that the diffusing walkers are extended objects — the PSPs in case P have 
finite size, and the domain walls in case Q have finite range of repulsion between them. In 
spite of this, in both cases the dynamics can be mapped exactly to the well known simple 
exclusion process (SEP) JTIJ of particles on a line of a reduced length <Lx = fxC- ■ Here £ 
differs slightly from £, and is defined below. 

In case P, for a background of type (a), every PSP is mapped onto a single particle, 
while a doublet is mapped onto a hole (TTIITTTTTI4TTTTT-U - * ° ° • ° *°) leading to r x = 
(1 — 3p)/2. In the case of backgrounds of type (b), which have only a single PSP, £x depends 



on the number of || pairs present; if the maximum number of such pairs are allowed (e.g. 
TTITI4T4TT •••)> r x = 3/4, as every doublet and a single spin gets mapped to a hole. In case 
Q, every parallel triplet is replaced by a hole and an f j pair associated with a domain wall is 
replaced by a particle (TTTTIIIIIIIITTTT - > °«°°«°) leading to r x = (l+p)/3. The modified 
length £ is (£ + 4) and (i + 2) in the cases P and Q respecively The difference between 
£ and £ arises from consideration of the dynamics of the walkers near the active- quiescent 
boundaries; since £ ~ £ for large £, we have p = p£/£ ~ p. 

For the mapped SEP with free boundary conditions, one knows that the relaxation time 
— ^-x /t^ 2 D where D = 1/2 is the diffusion constant of the SEP. Hence in our problem 
Tx = 2r x £ /tt 2 . The full autocorrelation function C x (t) is then given by averaging over the 
distribution A(£, p) of stretches of length £ and walker density p. 

C x (t)=J2A(£,p)C £ x (p,t) (2) 

e, P 

To make further progress, note that the largest and the smallest values of tx (as p is varied 
from l/£ to its maximum possible value) provide bounds on Cg, (p, t), and thereby on C x (t). 
Using A(£) ~ exp(—\x£) and r x = 2r* x £ /tt 2 , (where r x is a different constant for the 
minimum and maximum densities), we see that the above sum is dominated by a saddle point 
value £* = {-kH / \xr* x )^ ■ This leads to a stretched exponential form for both the upper and 
lower bounds, implying that C x (t) ~ exp(—{t/T X )^). In case P, the lower and upper bounds 
on To are (l/5) 2 (8/(277r 2 Ap)) (corresponding to p = 1/5) and 8/(277r 2 Ap) (corresponding to 
p = 0). In case Q, they are (l/2) 2 (8/(27vr 2 A 2 3 )) (corresponding to p = 1/2) and 8/(27vr 2 A|) 
(corresponding to p = 0). Inserting numerical values Ap = 0.22 and Aq = 0.13, we have 
0.024 < r^ < 0.61 and 0.44 < r® < 1.75. This helps us to understand the reason for the 
numerically observed slower decay (r^ = 0.37 and t® = 0.56 in Fig. 4) of C Q (t) compared 
to C p (t). 

4 Conclusion 

In our study ofT = oo— > T = quenches of the extended ANNNI model evolving under 
triplet flip dynamics, we found interesting arrested states with interspersed quiescent and 
active regions. Although the character of excitations is quite different in the two cases 
studied, the autocorrelation function shows similar stretched exponential decays. This was 
traced to the fact that active regions in both cases carry diffusing excitations, and the 
dynamics in both cases can be mapped on to the simple exclusion process. 

We have found and studied IQA states in quenches of other models including the ANNNI 
model with Kawasaki dynamics. In this case, in addition to the stretched exponential decay 
of C(t) in steady state, we found an interesting 2-step relaxation towards equilibrium when 
the temperature is raised slightly JTT| . 

It would be interesting to try to identify signatures of IQA states in rapid quenches of 
polytypic materials, for instance near the 3C — 6H phase boundary. 
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